Effects of imperfect noise correlations on decoherence-free subsystems: SU(2) 

diffusion model 
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We present a model of an JV-qubit channel where consecutive qubits experience correlated random 
rotations. Our model is an extension to the standard decoherence-free subsystems approach (DFS) 
which assumes that all the qubits experience the same disturbance. The variation of rotations acting 
on consecutive qubits is modeled as diffusion on the SU(2) group. The model may be applied to 
spins traveling in a varying magnetic field, or to photons passing through a fiber whose birefringence 
fluctuates over the time separation between photons. We derive an explicit formula describing the 
action of the channel on an arbitrary iV-qubit state. For TV = 3 we investigate the effects of 
diffusion on both classical and quantum capacity of the channel. We observe that nonorthogonal 
states are necessary to achieve the optimal classical capacity. Furthermore we find the threshold for 
the diffusion parameter above which coherent information of the channel vanishes. 

PACS numbers: 03.67.Hk, 03.65.Yz, 42.81.-i 



INTRODUCTION 



A fruitful approach to protect quantum systems from 
decoherence introduced by uncontrolled interactions with 
the environment is to use symmetries exhibited by those 
interactions. When elementary quantum systems in an 
ensemble are coupled to the environment in an identical 
way, it is possible to identify certain collective degrees of 
freedom that turn out to be completely decoupled from 
the interaction, thus preserving quantum coherence. This 
approach is the basic idea behind decoherence-free sub- 
spaces and subsystems (DFSs) 0, 0] (for a review see 
[3|), which can be implemented in a number of physical 
scenarios including quantum dots Q, atoms in a cavity 
[j| , ion traps @[ or photons transmitted through a bire- 
fringent fiber [3]. The existence of DFSs allows one to 
encode or communicate reliably both classical and quan- 
tum information, which has been demonstrated in first 
proof-of-principle experiments [H]. 

The purpose of this paper is to go beyond the stan- 
dard theory of DFSs and to analyze a scenario when 
the interactions of elementary systems with the environ- 
ment are not necessarily identical. Specifically, we will 
consider a sequence of qubits affected by random uni- 
tary transformations. In the standard DFSs model all 
unitaries are the same, enabling one to apply the stan- 
dard decomposition into multiplicity subspaces based on 
angular-momentum algebra. In this paper, we will as- 
sume that correlations between unitaries affecting con- 
secutive qubits are characterized by a probability distri- 
bution describing isotropic diffusion on the SU(2) group. 
Such a distribution has a number of invariant properties 
that result in the existence of a single real parameter 
which defines the strength of correlations. This will en- 
able us to investigate the continuous transition between 
the extreme regimes of perfect correlations and indepen- 
dent transformations of consecutive qubits. 



The model introduced in this paper is a certain channel 
acting on a sequence of N qubits. We derive its detailed 
description using angular momentum algebra, and ana- 
lyze quantitatively the case of N — 3 qubits. In particu- 
lar we present numerical results concerning the classical 
capacity and present the optimal ensemble of states at- 
taining it. Surprisingly, for imperfect correlations the op- 
timal ensemble contains nonorthogonal states. The fact 
that the use of nonorthogonal states may be necessary 
for the optimal classical communication has been known 
[y], yet it is rather surprising that this effect appears in a 
channel constructed in such a natural way. We compare 
the optimal capacity with the one that can be achieved 
when one is restricted to the use of orthogonal states 
and indicate orthogonal states which perform almost op- 
timal. We calculate and analyze the optimal coherent 
information, which provides certain information about 
the quantum capacity of the channel. 

Results presented in this paper can be applied to differ- 
ent physical systems used for quantum communication, 
e.g. spins traveling in fluctuating external magnetic field, 
and photons traveling in a fiber with randomly varying 
birefringence. The model generalizes the standard DFS 
theory whenever the time separation between qubits sent 
through the channel is not negligible in comparison with 
the characteristic time of environment fluctuations. One 
can easily imagine such a situation in the case of spins 
traveling relatively slowly in randomly varying magnetic 
field. When it comes to photons, however, in most of 
present experiments the standard DFS approach seems 
satisfactory. The time separation of photons is typically 
several orders of magnitude smaller than the fluctuation 
time of birefringence, whose variations are caused by me- 
chanical and thermal factors and hence by their nature 
are relatively slow. Nevertheless, our model might be 
relevant also in this setting, when large temporal sepa- 
rations between photons need to be introduced, or se- 
quences of photons are long enough to reach the time 
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scale of birefringence fluctuations. 

The paper is organized as follows. We describe the 
model of imperfect correlations in Sec. HH In Sec. ITTT1 
investigate the structure of output states and give an 
explicit and efficient formula for calculating the action of 
the channel on an arbitrary N qubit state. In Sec. |IV] 
we consider the example of three qubit communication, 
which is the smallest number of qubits allowing quantum 
information to be transmitted. Finally, Sec. IVl concludes 
the paper. 



where p(U) is a certain probability distribution defined 
on the group SU(2), which fully characterizes our deco- 
herence model. 

As the explicit model for the distribution p(U) we will 
take the solution to the diffusion equation on the group 
SU(2). This solution can be conveniently parameterized 
with a non-negative dimensionless time t characterizing 
the diffusion strength. The explicit expression for the 
distribution Pt{U) on the SU(2) group at time t reads 



II. IMPERFECT CORRELATIONS 

The standard model for collective depolarization of N 
qubits described by a joint density matrix p is given by 
the following map: 

T(p) = JdU U® N ptf® N , (1) 

where U is an SU(2) matrix describing the rotation of a 
Bloch vector of a single qubit state and dU is the Haar 
measure on the SU(2) group. This operation inflicts a 
rotation which is completely random and identical for all 
qubits. In other words we may say that noise experienced 
by a given qubit is perfectly correlated with the noise 
experienced by all other qubits. We will refer to this 
map as the twirling map. 

Let us now extend the above model in order to describe 
the situation where correlation of noise acting on different 
qubits is not perfect. We will assume that the density ma- 
trix p of the qubits undergoes a unitary transformation 
given by a tensor product U\®U2®- ■ -®Un averaged with 
a certain probability distribution p(U\, Ui, ■ ■ ■ , Un)- We 
will treat the sequence of Ui as a discrete-time stochastic 
process with the Markov property, i.e. as a Markov chain 
[nj]. Consequently the probability distribution for Ui will 
only depend on Ui-i. Additionally, we will impose the 
stationarity condition on our process which means that 
the conditional probability distribution is described by 
the same function irrespectively of the index i. Thus the 
joint probability distribution is given by a product 

p(U 1 ,U 2 , ...,U N )= p{U n \U N -i) ■ • 

(2) 

We will consider an isotropic process, i.e. one that does 
not distinguish any element of the SU{2) group, which 
is equivalent to the physical assumption that the fluc- 
tuations do not favor any particular form of depolariza- 
tion. The isotropy condition requires that the conditional 
probability distribution p{Ui\Ui-i) = p(VUi\VUi-i) = 
p{U i V\U i - 1 V), for every V € SU(2). This implies that 

KW-i) = KtWti) = p(uUUi), (3) 



Pi ([/)=£(2j + l)e x pf- 5 j(j + l)*) £ &(U)Z 

j—0 m——j 

(4) 

where Qi jU)™ are rotation matrices of the SU(2) group 
E3, EH- The function pt(U) is invariant with re- 
spect to unitary transformations of its argument, i.e. 
p t (VUV^) = p t {U) for any U, V G SU(2). It depends 
only on the eigenvalues of U which can be parameterized 
with a single angle £ as e 1 ^ 2 and e~^/ 2 . Then the sum 
over m in Eq. (J4j) , equal to the character of the jth irre- 
ducible representation, can be calculated explicitly as: 

V &(U) m - Sin[(j ± 1/m (5) 
L [ )m ~ sin(£/2) ' (b) 

m=—j 

The family of distributions pt(U) has the following con- 
volution property: 

J dU Ptl (U'U^p t2 (U)=p tl+t2 (U'). (6) 

For t — > the function pt{U) tends to a delta- like distri- 
bution peaked at identity, while for t — > 00 it becomes 
uniform on the group. Therefore the limit t — > cor- 
responds to perfectly correlated depolarization specified 
in Eq. ((TJ, and the opposite case of independent depo- 
larization affecting each one of qubits is recovered when 
t — * 00. We will be interested in perturbations in DFSs 
resulting from imperfect correlations characterized by a 
finite value of t, 



III. CHANNEL ACTION 

We now proceed to write explicitly the action of the 
channel on an arbitrary TV qubit state p. The reasoning 
below is valid for an arbitrary isotropic function p(U), 
and the exact form of p(U) specified in Eq. (J4j) is not 
used until Eq. lfT6|) . The output state of the channel 
reads: 
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■■■®U J N = 



£(p) = J dU! J dU 2 . . . J dU N p(Ux, U 2 ,..., U N )Ui ®U 2 ®---®U N P U\®UI 
dUt [ dU 2 ... [ dU N p{U 2 \U 1 )p{U s \U 2 ) . . .piUtfpN-jU! ® U 2 ® • • • U N p U\ g> u\ • • • ® C/^ 



Taking into account Eq. ([3]), and introducing new variables t/,' = Uillj_ 1 we arrive at: 



«ir, I du 2 ...J du' N p{u' 2 )p{u' a )...p(u' N )x 

U x <g> (U'zUi) <8 • • • (U'tfU'tf^ . . . U x ) p U\ (g> {u\u'^) ® • • • ® ([/Jk 



dC/^p(C/jv)(l <g> • • • ® £4) • 



t/T't 



d/7i(?7i ® 
I 



£7i) p (Ul <g> ■ ■ ■ ® Ul)) ...(1 



U%). (7) 



The action of the channel can thus be written in a com- 
pact form: 

e(p)=i N - 1 (i N - 2 (...i 1 (T( P ))...)), (8) 

where: 

Up) = J du P (u) 

x l(g> l(g>U (g> ■ ■ ■ (g>U p l(g>l<S)U^ ■ 

N-i i N-i 



U[. (9) 



It can be easily proven that thanks to the isotropy 
of the function p{U) specified in Eq. (J3j) , operations Xj 
commute with each other, i.e. Ii(Ik(p)) = 
and furthermore they commute with T. This implies 
that T(£(p)) = £{T{p)), which together with the fact 
that T(T(p)) = T(p) leads us to the conclusion that 
T(£(p)) = £(p). This means that the output state is 
invariant under the twirling map. We will say that it has 
a twirled structure. 

The action of the twirling map in Eq. ([8]) allows 
us to restrict our considerations to input states hav- 
ing a twirled structure. Let us now describe the states 
Ptwiried = T{p). For this purpose we use the following 
standard decomposition of the space of N qubits: 



n 



®N 



N/2 

e 

j — (N mod 2)/2 



Hi 



dj times 



N/2 



j=(N mod 2)/2 



There are different ways of obtaining a given angu- 
lar momentum by adding elementary 1/2-spins of sin- 
gle qubits. If a given angular momentum is obtained 
by always adding the z-th spin-1/2 to the previously ob- 
tained total angular momentum of i — 1 spins, these 
different ways correspond to different paths in the van 
Vleck diagram [TH], shown in Fig. [TJ Therefore a can 
be regarded as an index labelling different paths lead- 
ing to equivalent representation subspaces. There are, 
however, alternative ways of specifying a which depend 
on the order of adding elementary 1/2 spins. In the 
following calculations, it will be helpful to use differ- 
ent conventions for labelling equivalent representations, 
distinguished with a subscript {k}. In explicit terms, 
a{ fe } will label equivalent representation subspaces which 
are obtained by first adding successively the first k 
spins 1/2, from number 1 to k, then adding the last 
N — k spins in the reverse order, from number N to 
k + 1, and finally combining these two groups together. 
The value of a^k} is thus fully specified by a sequence 
{ji, 3i2, ■ ■ -,ji...k}{jk+i...N,jk+2...N, ■ ■ -,3n}, where j l ... k 
is the total angular momentum after adding together 
spins from the number i to k. Of course, a single par- 
ticle spin is equal to ji — 1/2 for any i. Naturally, all 
the conventions for numbering equivalent subspaces, cor- 
responding to a different choice of the k, are legitimate, 
yet the ability to switch between alternative conventions 
will make it possible to derive a compact formula for the 
-dj j action of the channel. 



(10) 

where TLj is a representation subspace in which the (2j + 
l)-dimensional irreducible representation of SU(2) acts, 
and Cdj is the multiplicity subspace which is not affected 
by the action of J7® . The above decomposition allows 
us to use a basis \j, m, a), where j is the total angular 
momentum, m is projection of angular momentum on 
the z axis, and a labels representation subspaces that 
are equivalent, i.e. have identical j. 



Using the decomposition given in Eq. lflO|) we can write 
any twirled state ptwiried in the following way: 



Ptwirled = T{p) 



N/2 

e 



j=(N mod 2)/2 



2.7 



jtn^Pj, (11) 
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FIG. 1: Van Vleck's branching diagram describing addition 
of N spins 1/2. The X and Y axis corresponds to number of 
spins N and the total spin j respectively. The numbers in cir- 
cles represent the multiplicity of given total angular momen- 
tum subspace (number of paths leading to it). An exemplary 
path {|, 1, i, 0, |, 1} is marked with bold (red) line. 

Denoting by Pj a projection onto Tij <£> C<j . we find the 
relations 

Pj = Tr(P iP ), (12) 
Pj = -'in, (PjpP 3 ). (13) 
I 



Because the action of Xj preserves the twirled structure, 
we may write the output state of the channel in an anal- 
ogous form: 

N/2 out 

£{ P ) = £(/w led ) = ^TT 1 ^ ® p 5 ut - 

j — (N mod 2)/2 

(14) 

Therefore, the full description of the action of the chan- 
nel amounts to writing p° ut ,p° ut as a function of Pj,pj- 
Equivalently we may characterize the channel by its ac- 
tion on operators: 

P J W = VT-T E \J,M,a){J,M,a'\ = 

ZJ + 1 m=-j (15) 

= -^-^Hj®\oc){a'\ 

since any twirled state can be written as their linear com- 
bination. As discussed before, we have freedom in choos- 
ing the convention for labels a. The simplest expres- 
sion is obtained if we specify how the channel acts on 

Pj , and express the output in terms of the opera- 

tors Pj {jv ~ 1} _ ihe detailed derivation is presented 

in Appendix El with the final formula in the form: 



£{P[ 



J 



31+32. ..N 



312+33. ..N 



31...N-1+3N 



£••■ E E 



E 



E 



p 



'{N-l},a {N _ 1} 
Jn-i 



312 
312 



31...N-1 J 1 = \31-32...n\ J2 = \jl2-33...N 
31...N-1 



n(4.\Ji,3i,32...N 



R(t) 



• / 2.Jl2.J3.. 
Jl,3l2,33.. 



Jn— l=\jl. ..N — l—jN 

J1...JV-1 J 
jN — l,N 

JiV-l,ji...jv-iJjtf 
JN — 2,jl...N — ldN 



u(j n . 2 ,n - i)^u(j N _ 2 , n - 1)^--; x 



(16) 



where according to our convention a>{k}i and a'| fc | read: 

«{fc} = {il)il2, ■ ■ ■ , jl...k}{jk+l,...,N, jk+2,...,N, ■ ■ ■ Jn} 
a '{k} = {il)Jl2) ■ ■ ■ i 3l...k}{3k+l,...,N> 3k+2,...,Ni ■ ■ ■ >3n}i 

(17) 

the functions R(t) j*^ '^' 2 are given in terms of Wigner 6j 
coefficients [l2] as: 

fl (*&^= E (-l) J ^ J (2j + l)(2J fe + l)x 

i=|j3— Jal 

(18) 



and f7(J, fc)^ '"^ is a shorthand notation for specific co- 
efficients t/(ji, J2j_/, J3; Ji2, J23) used in adding three an- 
gular momenta [l2|, discussed in Appendix [Aj* 

^ (-7. = U{j 1 ... k - 1 ,l/2,J,j k+ i... N ;ji...k,jk...N)- 

(19) 

Recall also that all single particle spins j'j = 1/2, which 
we kept implicit in Eq. lfl6|) only to ease the understand- 
ing of the formula structure. 

The formula (fT6|) . despite its lengthy appearance, en- 
ables an efficient calculation of the channel action even 
for a large number of qubits, which would be infeasible 
by inserting directly the formula for Pt(U) into Eq. J7]) 
and performing the integration. 
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IV. THREE QUBITS 

The case N = 1 of one qubit is trivial. The channel 
acts as completely depolarizing channel, hence its capac- 
ity both classical and quantum is zero. The case N = 2 of 
two qubits was solved in Ref. [16], where the optimal clas- 
sical capacity was derived together with the optimal en- 
semble of states. The twirling operation in this case pro- 
duces Werner states p = (1 + c)l/4 — c\ip~) |. Imper- 
fect correlations were introduced in that paper through a 
phenomenological shrinking factor 77 multiplying the pa- 
rameter c. Our general model predicts the same behav- 
ior, with the shrinking factor given explicitly by 77 = e - '. 
Notice also that quantum capacity of this channel is zero 
since all multiplicity subspaces are one dimensional. 

We will now investigate the action of the channel for 
three qubits, which is the lowest number with non-trivial 
equivalent representation subspaces. Adding three spins 
1/2 gives one subspace with j = 3/2, which is the fully 
symmetric subspace, and two subspaces with j = 1/2. 
The two-fold multiplicity of subspaces corresponding to 
j = 1/2 enables one to preserve quantum superpositions. 
In the case of perfect correlations t = the channel acts 
as the identity in the multiplicity subspace [2|, which 
therefore forms a decoherence-free subsystem, making it 
possible to encode one qubit with the fidelity equal to 
one. The classical capacity of such a channel is equal 
to log 2 3, since we have three completely distinguishable 
output states: two orthogonal states of a qubit from the 
multiplicity subspace and one state from the fully sym- 
metric subspace. When correlations between operations 
acting on consecutive qubits are not perfect the fidelity 
of the transmission as well as classical and quantum ca- 
pacity will decrease. In this section we will derive ana- 
lytical formulas for the fidelity of transmission through 
the channel as well as numerical and approximate ana- 
lytical results for quantum and classical capacities of the 

I 



where the output matrix is written in the basis \e{), |e 2 ), 
1 2). This expression will be now used to calculate quan- 
tities characterizing the channel. 



channel. 

For three qubits a general twirled state has the form: 

P=f(l« 1/2 ®Pi/2)©^--^ 3/2 , (20) 

where pi/ 2 is an arbitrary 2x2 density matrix, and < 
p < 1. Let us now specify a particular basis in the j = 1/2 
multiplicity subspace in which the operation £ will have 
the simplest form. Let the state |0) correspond to the 
equivalent subspace labeled {j = 1/2, a = 0}, obtained 
by combining the first two spins 1/2 together to get the 
angular momentum and then adding the third spin with 
the total angular momentum 1/2, while |1) corresponds 
to the equivalent subspace labeled {j = 1/2, a = 1}, 
obtained by combining the first two spins to get the an- 
gular momentum 1, and then adding the third spin to 
produce the total angular momentum 1/2. For simplicity 
in what follows we will omit the identity operators acting 
in the representation subspaces. Instead of tn 1/2 ® P1/2 
we will simply write a single qubit state P1/2, and instead 
of l« 3/2 we will write |2)(2|. Our channel can be thus 
regarded as effectively a qutrit channel that preserves no 
coherence between subspaces {|0),|1)} and {|2)}. The 
evolution of the state has the simplest form if we intro- 
duce the following basis in the qubit subspace: 

|0) + V3|1) V3|0)-|l) 
l e i) = 2 ' ' 2 ' = 2 ' ^ ' 

In order to describe the action of the channel it is suffi- 
cient to calculate its action on the extreme states whose 
convex combinations generate the entire set of twirled 
states. The extreme states are |2)(2| and an arbitrary 
pure state \ip)(i>\ in the qubit subspace, which we will 
parameterize as \ip) = cos(0/2)|ei) +sin(6>/2)e^|e 2 ). Ap- 
plying Eq. (fl6| to the three-qubit case yields: 



(22) 



(23) 

I 

A. Fidelity 



Let us first calculate the transmission fidelity of a pure 
qubit state encoded in the multiplicity subspace. Accord- 
ing to Eq. f23j). with a probability |[3- 4e~* sin 2 (6»/2) - 
e _2 *(l + 2cos#)] the input qubit state is removed from 



WD=( 1(l ei £(3-4e°-« + e-*) ) ® ^ + ^ 4 



e~ 2t ) 



£{m{ ^ ) _ l . l{l + e- 2t (l + 2co S 6)) isin0( e -*cos 



\ sin0(e-* cos + ie- 2t sin cj>) i (3 + 8e~* sin 2 (<9/2) - e - 2t (l + 2 cos 6)) 

© i (3-4e-*sin 2 ((9/2)-e- 2 *(l-|-2cos(9)) , 
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the multiplicity subspace and transformed into the state 
|2)(2|. We can consider an effective one qubit channel 
£efi by replacing the state |2)(2| at the output with the 
maximally mixed state in the qubit space, given by 1/2. 

Since now we deal with an ordinary one qubit channel, 
i.e. a one qubit trace-preserving completely positive map, 
we may investigate its action in terms of an affine map 
on the Bloch vector. Written in the basis |ei), fa), the 
output Bloch vector r out can be expressed in terms of the 
input vector r; n as: 














2e~ 




(24) 

The channel shrinks and translates the initial Bloch 
sphere. Notice that the shrinking is not isotropic and 
is weakest and strongest respectively in the x and y di- 
rections. The translation magnitude initially increases 
achieving the maximal value at t = ln(2) and then van- 
ishes asymptotically. 

The explicit expression for the fidelity of the output 
state is given by: 

f(6,<j>,t) = (VI |V) - |(1 + r in .r out ) = 

= — [12 + 5e-* + 7e~ 2t + 
24 

+ (e^ 2 * - e _t ) (4 cos 6 - 6 cos 20 sin 2 9 + cos 29)] . (25) 

and is depicted in Fig. [21 It is seen that the fidelity is op- 
timized along meridians cj> = 0, n, which is a consequence 
of the fact that shrinking is weakest in the x direction. 
T he h ighest t rans mission fidelity is achieved for states 
i/5/8|ei) ± y/Z/%fa), independently of the actual value 
of the diffusion time t > 0. 

We now compute the average fidelity of states lying on 
a great circle parameterized by the unit normal vector 
given in spherical coordinates (6 c ,4>c)- This quantity is 
relevant when one needs to transmit a relative phase of an 
equally weighted superposition [T3|. The average great 
circle fidelity reads: 



/c(0c,0c) = 

- 2t + e~* + 2) + (e- 2t - e -t ) (1 + 3 cos 2</> c ) sin 2 



-[6 
24 L 



(/(*)> 



1 
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Thus it is a straightforward observation that the great 
circle parameterized by (9 C — 7r/2, cj) c = ir/2) (the merid- 
ian in the xz plane) is optimal and yields f c = (3 + 2e~' + 
e- 2t )/6. 

The average fidelity integrated over the entire Bloch 
sphere of input states with the uniform distribution is 
given by: 



FIG. 2: The fidelity f{6,4>,t = 1) of states transmitted 
through an effective qubit channel £ e g as a function of the 
Bloch parameters 8 and <j> of the input state. 



and it decreases monotonically with increasing diffusion 
strength. In the limit of strong diffusion t 1 the the 
channel is completely depolarizing. 



B. Coherent information 

We will now optimize the coherent information of the 
channel £ in order to estimate its capacity for transmit- 
ting quantum information. Let us note that we now need 
to consider the full qutrit channel £ acting according to 
Eqs. (|22|) and (|23|) rather than the effective one qubit 
channel £ e ff, as considering only the latter would lower 
the achievable quantum capacity. 

Coherent information is defined as follows [18]: 

I c = sup (S(£(p)) - S env (£, p)) , (28) 
p 

where S(p) = — Tr (plog 2 p) is the von Neumann entropy 
and S'env is the entropy exchange [l8j]. The analytical 
optimization seems hard due to the complicated form of 
the Krauss operators [18] our channel has. Nevertheless, 
c] one can easily notice that the optimal state p in Eq. f28|) 
(26y m be supported on |Jei), fa)} subspace, since the sym- 
metric subspace \2) can be regarded as a purely classi- 
cal degree of freedom and therefore cannot contribute to 
quantum capacity. Numerical optimization shows that 
the optimal state has the form: 



p(e) = e|ei>(ei| + (l-e)|e 2 )<e 2 |. 



(29) 



(9 + 4e~ 



5e 



~2t\ 



(27) 



We have optimized analytically the coherent information 
in the limit of weak diffusion (small t). In the case of 
no diffusion the state p(e =1/2) maximizes the quantum 
capacity: Ic = 1. We have expanded Ic in a power series 
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FIG. 3: Coherent information Iq as a function of diffusion 
time t. The numerical solution is marked with the solid line, 
the analytical approximation given by Eq. (|30|> is marked with 
the dashed line. 




0.05 0.1 0.15 0.2 0.25 0.3 



FIG. 4: The parameter e of the input state maximizing the 
coherent information Ic- The numerical solution is marked 
with the solid line, the analytical approximation given in Eq. 
([31]) is marked with the dashed line. 



in e around e = 1/2 up to the second order and calculated 
the maximum. In the lowest order in t (which includes 
also terms tlogt, reflecting the fact that the derivative 
of Ic(t) diverges in t — 0) the approximate input state 
parameter e and the quantum capacity read: 



Ic 



1 - ^ ( • - . • 



In 2 



- 2 log 2 t 



log 2 3 



(30) 
(31) 



Numerical results depicted in Fig. [3] indicate that coher- 
ent 

information drops to strictly zero for diffusion time 
t > 0.275, which strongly suggest that in this regime no 
quantum communication is possible. Since our channel 
is not degradable [19( it could happen that the quantum 
capacity is not zero even though the coherent information 
vanishes, see e.g. Ref [20l |. 



C. Classical capacity 

The classical capacity of a quantum channel can be 
calculated using the Holevo-Schumacher- Westmoreland 
GH HI formula: 



C = 



sup 

{pi.Pi} 



S £ 



Pi Pi 



(32) 

where the supremum is taken over all ensembles {pi, Pi}. 

Let C = S {£ {YstPtPi)) - J2iPiS(£ (Pi)) denote the ex- 
pression that is optimized in Eq. (|32| . In order to achieve 
the supremum it is enough to use d 2 pure states, where 
d is the dimension of the input Hilbert space [22|. In 
our case the maximum number of pure stares needed to 
achieve the supremum is five: four in the qubit subspace 
and additionally the symmetric state |2). We imple- 
mented numerical optimization of the classical capacity, 
which for all diffusion times t yielded three-element opti- 
mal ensembles containing only two states from the qubit 
subspace. Furthermore, the sub-ensemble in the qubit 
subspace was composed of states given by: 



|V>i> =cos(0/2)| ei )+sin(0/2)|e 2 ) 
IV^) =cos(0/2)| ei >-Bin(0/2)|ea) 



(33) 



and equal weights, which we will denote by q. The pa- 
rameters q and 9 are functions of t. 

For t > the optimal states \tpi) and \ip2) turn out to 
be non-orthogonal. In order to judge how significant this 
nonorthogonality is for the channel transmission, let us 
compare the optimal capacity with the capacity attain- 
able when composing the input ensemble from a pair of 
orthogonal states in the qubit subspace and the symmet- 
ric state 1 2). The optimal orthogonal states for trans- 
mission in the qubit subspace are l/\/2(|ei) ± \e2)). As 
seen in Fig. [SI the capacity obtained using these states is 
very close to the absolute optimum. However, a different 
choice of orthogonal states can significantly deteriorate 
the capacity, with the worst case corresponding to the 
states 1/V2(\ei) ± i\e2)). The best and the worst per- 
formance of these pairs can be explained intuitively by 
the anisotropy of the effective qubit channel £ e g: the op- 
timal orthogonal states lay on the x axis characterized 
by the weakest shrinking, while the worst ones on the y 
axis where the shrinking is strongest. This argument is 
of course not rigorous, as we are dealing here with the 
full qutrit channel rather than the effective one £ e s in 
the qubit subspace. While in the limit t = of no diffu- 
sion any pair of orthogonal states performs equally well, 
this invariance is broken for t > owing to the asym- 
metric character of the channel. However, the particular 
choice of the orthogonal states given by 1/V2(\ei) ± \e2)) 
performs close to optimum for any diffusion strength. 

In order to obtain an approximate analytical solution 
we expand C = C(9,q,t) in a power series up to the 
second order in (8,q) around (tt/2, 1/3) and optimize an- 
alytically. This yields an approximate solution for a weak 
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optimal capacity 
best orthogonal states 
worst orthogonal states 
approximate formula 



0.05 0.1 0.15 



0.25 0.3 




0.8 
0.6 



FIG. 5: The classical capacity of the quantum channel C as a 
function of diffusion time t. The solid line depicts numerical 
results, the dashed line depicts the approximate analytical 
solution from Eq. (06]), while dashed-dotted and dotted lines 
represent capacities obtained using respectively the best and 
the worst pairs of orthogonal states. Notice that using proper 
orthogonal states allows one almost to achieve the optimal 
capacity. 



diffusion, given by: 

1 * 
q * 3 + 108 



7 



41og 2 3 



Int 



|-±(l-bg a 3 + 21nt) 



(34) 
(35) 



A comparison between numerical and approximate an- 
alytical solutions is shown in Fig. [6] The graphs show a 
good agreement of analytical solution with numerical re- 
sults. Substituting Eqs. (J34j) and ((35]) into (7(0, q, t) and 
retaining leading terms yields an approximate formula 
for the capacity: 

t ( 14+llln3 \ , . 

C « log 2 3 + - f 1 + 71og 2 tj , (36) 

which is compared with numerical results in Fig. [5j 



V. CONCLUSIONS AND DISCUSSION 

In this paper we have introduced a model an N qubit 
channel with imperfectly correlated noise, i.e. rotations 
inflicting consecutive qubits are subjected to the process 
of diffusion. We have given an explicit formula for the 
action of the channel on an arbitrary N qubit state and 
for TV = 3 we calculated the optimal classical and quan- 
tum capacities and the states which are optimal for com- 
munication. Interestingly, we have observed that clas- 
sical capacity is maximized when using nonorthogonal 
states. Additionally, we analyzed the robustness of dif- 
ferent orthogonal states, which in the case of no diffu- 
sion (perfect noise correlation) are optimal for classical 




OT05 0.1 0.15 0.2 0.25 0.3 



FIG. 6: The optimal weight q and the optimal Bloch param- 
eter 6, given in Eq. {33]), of input states belonging to the 
j = 1/2 multiplicity subspace maximizing classical informa- 
tion transmission. The solid lines depict numerical results, 
dashed lines represent approximate analytical solutions, given 
in Eqs. l134[| and l13"i[l . Notice that for t > the states are 
nonorthogonal. 



communication. We have indicated the most robust or- 
thogonal states which perform almost optimal for all dif- 
fusion times. We have also found a threshold for the 
diffusion time above which the coherent information is 
zero and hence most probably no quantum information 
can be transmitted. 

The model is a very natural extension of the stan- 
dard DFS theory, and it can be applied to any physi- 
cal situation when the action of the environment on the 
qubits can be described by a stationary Markov chain of 
SU(2) matrices, and where the transition probability is 
described by an isotropic diffusion process on the SU(2) 
group. This is equivalent to the assumption that consec- 
utive matrices are the result of an isotropic random walk 
on the SU(2) group (see [23] for a review on random walks 
on groups). The assumption can be justified for systems 
we discussed in the introduction: spins traveling in the 
presence of randomly varying magnetic field and photons 
transmitted through a fiber with birefringence fluctua- 
tions. In both cases consecutive qubits (spins or pho- 
tons) experience varying rotations in the channel, caused 
by fluctuations of the magnetic field or by birefringence. 
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These multiple small random contributions lead to the 
random walk of the effective matrix describing the ac- 
tion of the channel on consecutive qubits. If these fluc- 
tuations are isotropic then obviously the random walk is 
isotropic and it leads to the model we have presented. 
Furthermore, even if fluctuations are anisotropic the ef- 
fective random walk will be isotropic thanks to the as- 
sumption that the channel is long enough to completely 
depolarize every single qubit transmitted. This remark 
is relevant for fibers where the birefringence fluctuations 
are typically modeled as anisotropic [24( , with the princi- 
pal axis corresponding to a linear polarization. However, 
if such a fluctuation occurs at some intermediate point 
of a sufficiently long fiber, random polarization rotations 
introduced by the preceding and the following sections of 
the fiber should make such fluctuations isotropic. 
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APPENDIX A: DERIVATION OF THE ACTION 
OF THE CHANNEL £ ON OPERATORS Pj' a ' 

According to Eq. ((Sj) the action of the channel can be 
expressed as a composition of operations Ji defined in 
Eq. © . Let us calculate the action of li on operators 
P"' a introduced in Eq. (fT5l) . Since li acts non trivially 
only on the last N — i qubits it will be convenient to num- 
ber equivalent representation subspaces using the con- 
vention described in Sec. IIIH and write explicitly — 
{ji, 3i2, ■ ■ - ,ji...i}{ji+i,...,N,ji+2,...,N, ■ ■ -Jn}- The first 



step is to calculate the action of 2j on -Pj W M : 



t®t®u 



dU Pt (U)x 



U P 7 



t <g> (g) • • • <g) erf . 



N-i 



N-i 



(Al) 



Using a.^} for labelling equivalent representations allows 
us to decompose the state \J,M, ct{i}) using Clebsch- 
Gordan coefficients, denoted here with square brackets, 



according to: 



31+1.. .N 



Jl... 

m l. 



3i+l...N J 

m i+ i... N M 



\J,M,a {i} ) = 

mi...i=— ji. 

mi.. a, ai) ® \ji+i...N,m i+ i,,, N ,a 2 ) 



i first qubits 



N—i last qubits 



(A2) 



The total angular momenta of the first i and the last 
N — i qubits, given respectively by ji...i and ji+\...N are 
uniquely determined by a^. Similarly, the labelling of 
equivalent subspaces a\, a 2 within each block of qubits 
is determined by via the following relations: ol\ = 

{jl, jl2, • ■ ■ , a 2 = {ji+l}{ji+2...N, ■ ■ ■ Jn} 

(remember that ji = 1/2 for any i). An analogous de- 
composition can be applied to |J, M, a'uy). 



The next step is to write the operator P 



ing the decomposition given in Eq. I|A2|1 . The operation 
li acts trivially on the first i qubits, while the action of 
jj®(n-i) Qn last N — i qubits can be written with 
the help of Wigner rotation matrix 33- ? (f/)™ . By sub- 
stituting the explicit form of pt{U) given in Eq. |(4]) into 
Eq. (|Aip and using properties of Clebsch-Gordan coef- 
ficients, Wigner rotation matrices, and 6j Wigner sym- 
bols (see chapters 3,5,7 in Ref. [13]) we arrive after some 
lengthy calculations at a compact formula: 



TAPr 



where: 



31...i+3i+l...N 

E 

A=\h...i—ji+i...i 



tt \ l 'Jdl...i,3i+l...N ^Ji ' 

(A3) 



32+32 

^ (-l) J '- / (2j + l)(2J l + l)x 

J=b'a— Jal 



exp 



-Jti + 1)* 



jl J ji \ f jl j'l j 

32 3 j'l J 1 j' 2 h Ji 



(A4) 



and the curly brackets denote Wigner 6j symbols. Notice 
that the expression derived in Eq. l|A3j) can be used only 
if the index i in the operation li is identical with the in- 
dex i specifying the convention for numbering equivalent 
subspaces a^y If we want to apply Eq. (j A3|l to calculate 
the full action of the channel, given by 



£(p? a )=i N - l (...MHPj' a ))- 



it is convenient to start from -Pj {1> ' {1> for the represen- 
tation of the input state, and then to adjust the conven- 
tion after each step. To carry on this procedure we need 



(A5) 



to be able to express operators P 



«{i-i}.«{i-i> . 



in terms of 



P 

The necessary expression can be derived using the 
standard formalism for adding three angular momenta. 
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Consider three spins j±, 32, and j'3. One can write a 
basis using the total angular momentum of all spins in 
two different ways depending on the order in which spins 
were added together. A ket \ji, {3233)323, J, M) corre- 
sponds to a state with the total angular momentum J 
and projection on the z axis M, when spins 32, ja were 
first coupled together yielding the angular momentum j'23 
and finally the spin ji was added resulting in the total 
angular momentum J. Analogously {(3132)312,33, J, M) 
corresponds to the situation when first spins ji and ji 
are added and subsequently the spin j 3 joins them. A 
unitary operation U(j\, 32, J, 33; J12, J23) that transforms 
between these bases: 

ji+h 

\ji, {3233)323, J,M) = ^2 u (h,h, J, 33; 312, 323) 

ji2=\h—ja\ 

x\{jih)3i2,33,J,M) (A6) 
can be expressed using 6j Wigner symbols as: 



where the coefficients U (J, fc)^ are defined in Eq. (fL9|) . 
For completeness let us also write an inverse relation al- 
lowing for the lowering of the index i; 



p?««« = j2u{ J ,i)t^Pj {i - 1} '^u{j > 4- 



3i.,.N 
3i...N 



(A9) 



Equipped with the above formulas we can now calcu- 
late the action of the complete channel. For the input 

state expressed as a combination of operators Pj 111 ' {1> , 
the action of the operation X\ is, according to Eq. (|A3[) 
given by: 



, J1+J2...JV , 



Ji=|ii— J2...n\ 



(A10) 



U {31,32, J , 33; 312, 323) = v(2ii2 + l)(2j23 + 1) 

x r^-ih+h+J+h) { 3i 32 312 1 (A?) 

1 33 J 323 



In order to calculate the action of I2 we need to repre- 



sent P 



Eq. (JA8J) yields: 



in terms of operators Pj 



«{2}."{2} 



Using 



Specializing these general formulas to our calculation 
(i.e. replacing ji with ji...i-i, 32 with 1/2 and j'3 with 
ji+i,...jv) we can write: 



3l...i 
31.. A 



(A8) 



p 



->'^> = J2u{J,2)t, N P^' a - 

j'll 
312 



*U{J,2) J r . 



(All) 



We may now apply the operation Z2, whose action on 

the operators Pj i21 ' f2} is again given by Eq. (|A3p . The 
combined action of I2 and T\ thus reads: 



i 2 (ii(p; {i> " {i, )) = 



31+32...N 

E 

Ji = \ji-32...n\ 



t>(+\ Ji ,3\<32...N 



i'12 



U{Jl:^_. N U{Jl,2)^ 



312+33.. .N 

E 

J2 = \jl2-33...N\ 



R(t) 



J2,fl2>33. 



P 



«{2},a { 2} 



Jl,312,j3...N J2 



(A12) 



Iterating this procedure yields the explicit formula for 
the action of the channel £ given in Eq. (fTo) . Notice 
in this formula the output state is expressed in terms 
of operators p" {N " 1} ' a i N -^ _ jf desired this expression 



can be converted back to the representation of operators 
p«{i}< Q {i} k y applying repeatedly Eq. I|A9|) . 
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